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We study Borromean 2n-halo nuclei using effective field theory. We compute the universal scaling 
function that relates the mean-square matter radius of the 2 n halo to dimensionless ratios of two- 
and three-body energies. We use the experimental value of the rms matter radius of 22 C measured by 
Tanaka et al. (2010) [3] to put constraints on its 2 n separation energy and the 20 C — n virtual energy. 

We also explore the consequences of these constraints for the existence of excited Efimov states in 
this nucleus. We find that, for 22 C to have an rms matter radius within 1 — a of the experimental 
value, the two-neutron separation energy of 22 C needs to be below 100 keV. Consequently, this 
three-body halo system can have an excited Efimov state only if the 20 C — n system has a resonance 
within 1 keV of the scattering threshold. 
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In the last twenty years several nuclei where the neu- 
tron distribution extends far beyond that of the protons 
have been discovered. These nuclei, such as 11 Li and 
12 Be thus have neutron “halos” (see Refs. m for early 
reviews). Understanding the way in which nuclear struc- 
ture changes in these neutron-rich systems may provide 
important clues to the behavior of nuclei far from the 
N = Z line. 

The most neutron-rich isotope of Carbon yet produced, 
22 C, has recently been identified as another example of 
a halo system. Tanaka et al. measured the reaction 
cross-section of 22 C on a hydrogen target and, using 
Glauber calculations, deduced a 22 C rms matter radius of 
5.4 ±0.9 fm [3] . Their measurement implies that the two 
valence neutrons in 22 C preferentially occupy the lsi /2 
orbital and are weakly bound [3 . This conclusion is sup- 
ported by data on high-energy two-neutron removal from 
22 C |4]. Since 21 C is unbound [5], this suggests that 22 C 
is an s-wave Borromean halo nucleus with two neutrons 
orbiting a 20 C core. 

In this work, we implement such a description of 22 C 
in an effective field theory (EFT) for systems with short- 
range interactions. This EFT was developed to study 
few-nucleon systems with scattering length much larger 
than the interaction range (see Refs. jam for reviews). 
It was then applied to the a — n system in Refs. [Si] , ex- 
tended to study three-body systems in Ref. m , and ap- 
plied to 2n-halo nuclei in Refs. HUE]. At leading order 
(LO) in the theory the only inputs to the equations that 
describe a 2 n halo are the energies of the neutron-core 
resonance/bound-state, E nc , and the nn virtual bound 
state, E nn , as well as the binding energy, B , of the halo 
nucleus * 1 . All other properties of the nucleus are pre- 
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1 The binding energy of 22 C treated as a three-body system is 

equal to the two-neutron separation energy of the nucleus, in- 


dicted once these (together with the core-neutron mass 
ratio A) are specified. Here we use the EFT to com- 
pute the function that describes the mean-square matter 
radius of an arbitrary Borromean 2 n halo. 

We then focus on exploring properties of 22 C, treated 
as a 2 n- halo nucleus with a 20 C core. For 22 C, neither 
E nc nor B is well-known mm- We therefore use EFT 
to find constraints in the ( B , E nc ) plane using Tanaka et 
al.’s value of the rms matter radius. A similar strategy 
was recently pursued in Ref. m, although there a sim- 
pler model of 22 C’s structure was employed. The connec- 
tion between the binding energy and several low-energy 
properties of 22 C, including the rms matter radius, has 
also been explored in a three-body model by Ershov et 

al. [15] . 

Yamashita and collaborators investigated such correla- 
tions in halo nuclei already in 2004 m (see also Ref. m 
for a review). In 2011 Yamashita et al. [21] attempted to 
apply EFT to analyze the experiment of Ref. [3] . How- 
ever, as we shall discuss further below, an additional as- 
sumption was made in Ref. [21, which renders the re- 
sults of Yamashita et al. model-dependent. The results 
we obtain here therefore differ — both in principle and in 
practice — from those of Ref. [21] . 

The EFT description of three-body systems developed 
in Ref. [10] also provides insights into the Efimov ef- 
fect [22 (see Ref. [23] for a review). This phenomenon oc- 
curs in the three-body system when the two-body inter- 
action generates a scattering length, a, that is large com- 
pared to its range. For sufficiently large a there is a se- 
quence of three-body states whose properties are related 
by discrete scale transformations. Refs. HU CE2 M M 
explore the possibility of finding excited Efimov states 
in several 2n-halo nuclei, including 20 C and 22 C. We use 
our constraint on the binding energy of 22 C to discuss 
the possibility that such states occur there. 
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Effective field theory is a powerful tool to study few- 
body systems at low energies, because long-distance 
properties of the system are insensitive to the details 
of the underlying short-range interactions. In the case 
of the effective field theory for halo nuclei ( “Halo EFT” ) 
the breakdown scale of the EFT, A 0 , is of the order of the 
inverse of the range of the two-body interactions, and this 
is much larger than Q, the generic low- momentum scale 
of the system. Q can represent either the momentum of 
the process p, or 1/a. At LO, we ignore terms suppressed 
by Q/A 0 , which amounts to approximating the two-body 
potentials to be zero-range. The results of LO EFT calcu- 
lations are, therefore, universal in the sense that they are 
independent of the short-distance physics. The “range 
effects” can be systematically taken into account order- 
by-order in the Q / Ao expansion [251433] . However, even 
in an LO calculation, one can estimate the uncertainty 
due to these neglected higher-order terms. This feature of 
EFT allows us to account for its theoretical uncertainty, 
and makes it uniquely suited to constraining B and E nc 
of 22 C from the matter-radius measurement. 


FADDEEV EQUATIONS 


The two-body virtual energies, E nx , are related to the 
scattering lengths, a nx , by 


E — 
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where, x = n or c, fi nx is the reduced mass of the cor- 
responding two-body system, and the ellipses indicate 
higher-order corrections. Throughout this Letter, we 
work in units with ft = c = 1. 

Following Refs. uni mi, we write the Faddeev equa- 
tions for the spectator functions, F x (q ), which describe 
the relative motion of the spectator particle, x, and the 
center of mass of the other two particles for an s-wave 
three-body bound state of two neutrons of mass m each, 
and a core of mass Am. 


F n (q) = \ J q 3 /2 J d (M) 

{[Go ( Tr(q f ,q),q'-,B)+H(A )] t n {q ' ; B)F n {q') 

+G C o (n 2 (q,q , ),q'-,B)t c (q , ;B)F c (q')}, (2) 
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F c (q) = j dq' q ’ 2 j d{q.q') 

GS Mr, q),q'; B) t n (q B)F n (q'), (3) 

where the momenta 7r, dq, and dq are given by 
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The three-body Green’s functions G\\ and Gq are given 
by 
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The two-body t-matrices are 
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To solve the Faddeev equations, an ultraviolet cutoff, 
A, needs to be introduced to the integrals. The three- 
body contact interaction, H( A), is then required to can- 
cel the cutoff-dependence [TCt. H{ A) scales as mh( A)/A 2 
where ft, (A) ^ 0(1). We have verified that one three- 
body counter term, added to Gq in Eq. Q, is sufficient 
to renormalize Eqs. © and ©>• Furthermore, all the 
results in this Letter can also be replicated by adding a 
three-body counter term to each of the Green’s functions 
in Eqs. © and ©• 

In Ref. DU, Gaussian regulators were introduced in 
the two-body potential to regularize the effects of the 
short-range interactions on the two and three-body ob- 
servables at a cutoff scale A. This is equivalent to our 
zero-range approximation up to corrections of order Q / A 
in the EFT. In Refs. Fl9ll2Tj a subtraction was performed 
on the equations for the spectator functions. The sub- 
traction functions were then set by assuming that the 
Born approximation holds for an appropriate subtrac- 
tion point. This assumption is, however, not valid for 
the integral equations which are employed here USES]- 

The three-body wave function in the Jacobi represen- 
tation with the neutron as the spectator is 


*n(p,«) = Gg(p, q; B) {t n (q; B)F n (q) 

1 T l 

+ 2 J ^ d{p.q)t n (n' nn {p,q)-,B)F n (Tr' nn {p,q )) 

+tc(Kc(P,^B)F c (n! nc (p,q))}, ( 11 ) 
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and with the core as the spectator is 
*c(p,q)=G c 0 (p,q;B) {t c (q; B)F c (q) 

+ J d{p.q)t n (n' cn {p,q)-, B)F n (n' cn (p,q))}, (12) 

where the momenta 7?^ n , 7r^ c , and 7r' cn are given by 


We calculated mB(r q) at the unitary limit E nn = 
E nc = 0- It is plotted in Figure [l] as a function of A. 
These results are in qualitative agreement with those of 
Ref. [19], but our value of / is lower, which means that 
Yamashita et al. overpredict the unitary limit (r g) for 
a given B. The discrepancy is approximately 15% for 
A = 20. 


*nn(P> ?) = P ~ ((4Vt) & ( 13 ) 

q)=p - £ ( 14 ) 

and 

Tt'cn&q) =p- \q- ( 15 ) 

The one-body form factors are then given by 

B x (k 2 ) = j d pp 2 J dq q 2 J d (g • k^j 

d> x (p,q) 'Mp, \q- k\), ( 16 ) 

where the wave functions ^f x are normalized so that 
J~ x (0) = 1. The mean-square distance of the neutron 
from the center of mass of the core and the other neu- 
tron, (r 2 _ nc ), can be extracted from 

F n (k 2 ) = l- l -k 2 {r 2 n _ nc ) + ..., (17) 

and the mean-square distance of the core from the center 
of mass of the neutrons, (r 2 _ nn ), from 

B c {k 2 ) = 1 - ^k 2 (r 2 c _ nn ) + . . . . (18) 


A TWO-NEUTRON HALO NUCLEUS WITH A 
POINT-LIKE CORE 


Based on geometrical arguments, we obtain the fol- 
lowing formula for the mean-square matter radius of a 
two-neutron halo in the point-like core approximation, 

(ro 2 ): 


<^> = 


2(A+ !) 2 2 4A — 2 

(. A + 2) 3 v n ~ nc ' (A + 2) 3 V c_nn/ 


(19) 


At LO in Halo EFT, the quantity mB(r q) depends on all 
the variables featuring in the Faddeev equations: E nn , 
E nc: B and A. But, being dimensionless itself, it can 
only depend on dimensionless ratios of these four pa- 
rameters. Thus it is convenient to define the function 
/ ( E nn / F>, E nc / B ; A), as |19|: 


(20) 


0.20 

\ 

\ 

\ 



< 0.15 

\ 

\ 

\ 

\ 

\ 


: 

O 

0" 0.10 

\ 

V 

s 

\ 


: 





0.05 

- 


— 

0.00 





0 5 

10 15 

20 



A 


FIG. 1: 

The dimensionless 

function /(0,0;A), 

defined by 

Eq. (20), versus A. 




IMPLICATIONS FOR 22 C 

The function / can be calculated for any value of A, 
but specializing now to the case of interest for 22 C, Fig. [ 2 ] 
shows a three-dimensional plot of / ( E nn /B , E nc /B ; 20) 
in the ( E nn / , E nc / B) plane. The disagreement with the 
results of Ref. [21] appears to be worse at finite values of 
E nc and E nn than the 15% we found in the limit E nn = 
Enc = 0 . 


f(£ nn /B,£ nc /B; 20) 



FIG. 2: f (E nn /B,E nc /B] 20) versus (E nn /B, E nc /B). 


pig. m gives us the results we need in order to set a 
model-independent constraint on the binding energies of 
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21 C and 22 C. First, though, we must account for the fact 
that, when applied to this system, Halo EFT only reliably 
predicts the difference between the 22 C matter radius and 
that of 20 C. We account for the finite spatial extent of 
the core by including that effect in our expression for the 
mean-square matter radius of the two-neutron halo: 

( r2) = (r ° 2> + ZT2 (r2>core - (21) 

Here (r 2 ) core is the mean-square radius of the core, which 
we take from the 20 C rms radius of (2.98 =b 0.05) fm 
measured by Ozawa et al. [36]. In subsequent cal- 
culations we also use the value of E nn obtained from 
dun = (-18.7 ±0.6) fm [37]. 

To calculate the cutoff of our EFT, we approximate the 
range of the neutron-core interaction by the size of the 
20 C rms radius. We then estimate the relative error of our 
calculation by \JmE nn / Ao, \/2mE nc / A q or V2mB/ Ao, 
whichever is the largest. The spectrum of 20 C has also 
been measured [38]. It contains one bound 2 + state which 
lies 1.588 MeV above the ground state. We have not used 
this energy scale in assessing the breakdown of Halo EFT 
for 22 C, since the 2 + state can affect s-wave scattering 
processes only via higher-dimensional operators which do 
not enter the calculation at next-to-leading order. 

More generally, one might be concerned about the im- 
pact of neutrons in d-wave states on the structure of 22 C. 
The LO EFT being used here does not preclude the exis- 
tence of such states in either 20 C or 22 C, it only assumes 
that their primary effect on long-distance dynamics can 
be subsumed into the neutron-core and neutron-neutron- 
core contact interactions which appear in the leading- 
order calculation. 

Constraints on E nc and B 

In Figure [3} we plot the sets of (F>, E nc ) values that 
give v < r 2 > = 4.5 fm, 5.4 fm and 6.3 fm, along with 
the theoretical error bands. All sets of B and E nc values 
in the plotted region that lie within the area bounded 
by the edges of these bands give an rms matter radius 
within the combined (1-cr) experimental and theoretical 
error of the value extracted by Tanaka et ah. The figure 
shows that, regardless of the value of the 20 C — n virtual 
energy, Tanaka et al.’s experimental result puts a model- 
independent upper limit of 100 keV on the 2 n separation 
energy of 22 C. 

Since Yamashita et al. obtain an LO matter radius 
that is too large for a given binding energy, their con- 
straints on the maximum possible value of B are about 
20% weaker than ours. We also reiterate that their re- 
sults are based on an incorrect assumption regarding the 
high-energy behavior of the integral-equation kernel. Our 
result for the maximum binding energy of 22 C is a factor 
of two smaller than that found by Fortune and Sherr DU- 
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FIG. 3: Plots of (r 2 ) = 5.4 fm (blue, dashed), 6.3 fm (red, 
solid), and 4.5 fm (green, dotted), with their theoretical error 
bands, in the (H, E nc ) plane. 


Ref. [17] postulated that a simple extension of the model- 
independent relationship between B and (r^) that pre- 
vails in a one-neutron halo |39j applies to two-neutron 
halos. Our calculation suggests that the relationship pro- 
posed in Ref. PH does not accurately capture the three- 
body dynamics of the 2n-halo system. 


Implications for existence of excited Efimov states 

Following Refs. mum, we construct a region in the 
(H, E nc ) plane within which an excited Efimov state in 
22 C could occur. In Fig. [I]the purple region is that which 
allows at least one excited Efimov state above the ground 
state in an A= 20 Borromean nucleus. In the same plot, 
we also show the boundary curves that enclose the sets 
of E nc and B values which are consistent with an rms 
matter radius of 5.4 ± 0.9 fm once the theoretical errors 
are taken into account, i.e. those already displayed in 

Fig- m 

The Efimov-excited-state-allowed and rms-radius- 
constraint regions do not overlap for a 20 C — n virtual 
energy larger than a keV. (The 22 C radius can be com- 
puted accurately for 20 C —n virtual-state energies very 
close to threshold, but the computation of the existence 
of an Efimov state becomes numerically delicate here.) 
Indeed, as long as the trend in Fig. [4] continues, Efimov 
states seem to be precluded for values of E nc well below 1 
keV. However, we cannot make a stronger statement than 
this, since an Efimov state is present if we take E nc = 0, 
as a consequence of the fact that there are two 20 C— 77- 
pairs in the 22 C system [23]. Therefore, while we cannot 
categorically rule out the existence of an Efimov state in 
22 C, we can say that the 21 C system would need to be 
tuned very close to the unitary limit in order for one to 
be present. 
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FIG. 4: The region in the ( B , E nc ) plane that allows excited 
Efimov states (purple), and the region that encloses values 
consistent with the experimental rms matter radius of 5.4 zb 
0.9 fm (with same color-coding as in Figure [3]). 


CONCLUSION 


We used Halo EFT at leading order to examine the 
behavior of the rms matter radius of s-wave Borromean 
halo nuclei in the (E nn / B , E nc j B) plane. We computed 
the universal function / which describes this behavior in 
a model-independent fashion. We then applied these re- 
sults to 22 C, and put constraints on the (T>, E nc ) param- 
eter space using the experimental value of the 22 C matter 
radius. In contrast to previous works Hamm] which 
examined this problem, our constraint makes very few 
assumptions about the structure of 22 C. We use only the 
experimentally well-supported idea that it can be treated 
as a three-body system composed of 20 C and two neu- 
trons. Furthermore, our result incorporates the antici- 
pated theoretical uncertainty of the leading-order Halo 
EFT calculation based on this cluster picture. Even af- 
ter this uncertainty, and the experimental (l-cr) error bar, 
are taken into account we find that B < 100 keV for all 
values of E nc . This rules out the possibility of an excited 
Efimov state in the 22 C nucleus unless the 20 C — n system 
has a virtual state with an energy much smaller than 1 
keV. 
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